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WE will consider in the present paper a system of n differential equa- 
tions of the form: 


dy; = Pa, 7 
(1) 2 (te +4,;) 


in which the y,,’s are constants, and the a, ,’s are functions, not neces- 
sarily analytic, of the real independent variable z, continuous in the 
interval : 


We shall require that |a, ,| be integrable up to the point x = 0. 
For the development of certain sets of solutions we shall make the 
additional requirement that even after certain powers of log x have been 
multiplied into | a, ,| the resulting function shall be integrable up to the 
point z = 0; this requirement will be stated more explicitly later.} 

The point 2 = 0 satisfying these conditions may be called a regular 
singular point of the system of equations (1) in conformity with the use 
of that term by Professor Bécher in the study of linear differential 


* This paper was accepted in June, 1902, by the Faculty of Arts and Sciences 
of Harvard University in fulfilment of the requirement of a thesis fot the degree 


of Doctor of Philosophy. 
t The requirement that the functions a, , should be continuous in 0 < x < dis 


made only for the sake of simplicity. We might allow them to have a finite 
number of discontinuities in 06 of such a kind that each function |a, ,| can be 
integrated throughout the interval; and all the following work would hold with 
very little change. 

t Cf. p. 367. 
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equations of the second order.* ‘That this terminology is a legitimate 
extension of that commonly used when the coefficients of the system of 
differential equations are analytic functions of a complex variable, will be 
evident if the results of the present paper are compared with the thesis 
by Sauvage: Théorie générale des systémes d’équations différentielles 
linéaires et homogénes.t 

Our object is to investigate the nature of the solutions of (1) in the 
neighborhood of the regular singular point z = 0; and for this purpose 
we shall first reduce the system of equations to a canonical form by 
means of a linear transformation with constant coefficients of the de- 
pendent variable. We shall then apply the method of successive 
approximations to develop about the point x = 0 a system of n linearly 
independent solutions of the canonical system. By means of the linear 
transformation we shall return to ” linearly independent solutions of the 
original system ; and finally an application will be made to the case of 
the single homogeneous linear differential equation of the nth order. 


§ 1. 


A Spectat System or Equations: Its REDUCTION TO A . 
CANONICAL Form, anp So.ution.} 


Let us first examine the special case of (1) in which the coefficients 
a,, are all zero. In this case we have the system of differential 


equations : 


jon 
(3) 


where the yp, 7s are constants. 
A solution of this system may be obtained in the following way. 


Substitute 
= C; = constant 


* Cf. Trans. Am. Math. Soc., Vol. I. Jan. 1900, p. 41. The results of this paper 
are included as a special case in those we now give. Cf. § 7. 

+ Paris, 1895. Reprinted from the Annales de la Faculté des Sciences de 
Toulouse, Vols. VIIL. and IX. 

t The results of this section are not new, being on the one hand only slightly 
modified forms of Weierstrass’s results (cf. the foot-note on p. 345), and on the 
other hand special cases of the results obtained by Sauvage (cf. the last foot-note.) 
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in the equations (3), and then determine the constants r and C, so that 
the equations are satisfied. We obtain in this way the following system 
of n homogeneous linear equations for the C’s: 


(4) Bir Cr + O; (¢= 1, 2,...#). 


The necessary and sufficient condition that (4) may be satisfied by a 
set of C’s not all zero is that the determinant: 


Pie 


A(r)= 


. Pun 


shall vanish. This determinant equated to zero gives an equation of the 
nth degree in r, which is called the characteristic equation of (3);* the 
determinant itself we may call the characteristic determinant of (3). 

If the characteristic equation has distinct roots + 7,» We can 
determine 7 linearly independent solutions of (3) of the form : 


(6) = 0,527 1,2 


If however there is a multiple root, there will be, in general, solutions 
involving powers of log x. To determine these solutions, we must 
examine the minors of the characteristic determinant (5), and ascertain 
if this multiple root is also a root of all the first minors, second minors, 
etc. For the further study of this case, it will be useful to introduce the 
conception of tle elementary divisors of A(r). 

Suppose 7’ is a root of the determinant A (r) such that all the pth minors 
of A(r) are divisible by (r — rly, but no higher power of r — 1’ divides 
them all. In the same way (r — pl jord shall be the highest power of 
r — 1 dividing all the (p + 1)th minors. Then the expression : 


* Sauvage: c. p. 80. 
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is called by Weierstrass an elementary divisor of the determinant A (r).* 
It will be convenient to employ a different notation from that used in 
the definition of an elementary divisor. An elementary divisor of A (r) 
will be written : 


(r —_ T,) 
and it is to be noticed that several r’s with different subscripts, may be 


equal, as will be the case when a multiple root furnishes several elemen- 
‘tary divisors. We shall always have: 


k=m k=m 
(7) A(r) =] 
— | 


It can be shown (cf. the next foot-note) that a necessary and sufficient 
condition that a pair of systems of differential equations : 


j=n 
dx fit dx ja = 


can be transformed the one into the other by means of a transformation : 


> 45 2; 
j=l 


in which the A’s are constants whose determinant is not zero, is that the 
characteristic determinants of the two systems have the same elementary 
divisors. This theorem enables us to simplify the solution of the system 
(3); for we can write down a second system of differential equations 
having the same elementary divisors as (3), as follows: 


dz r, 


where 2. = 0. 


* Cf. Muth: Theorie und Anwendung der Elementartheiler, p. 2. 
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The characteristic determinant of (8) is: 


1 
1 
(9) A(r)= eves 
1 m—r 
1 +r,—r 


1 


and it will be easily seen that it has the elementary divisors (r — r,)%, 
(r ... Then, by the theorem above referred to, 
there exists a set of n? constants A,,,, whose determinant is not zero, 
such that: 

k=m 


(10) ¥: => (¢=1,2,...m). 


k=1 


The system of differential equations (8) we may speak of as the 
canonical system ; and now it is easily seen that: 

The canonical system of equations admits e, solutions, corresponding to 
the elementary divisor (r — 1,)°*, of the following simple form : 


k+k 
= l 

(log x) ASIS& 


and the n solutions obtained by giving « the values 1, 2,... m are linearly 
independent. 


* Cf. Weierstrass, Werke, Vol. II. pp. 75, 76. The case considered by Weier- 
strass is very easily reduced to the one we are considering by the change of inde- 
pendent variable t = log x. This reduction of Weierstrass is also given in Muth’s 
Elementartheiler, pp. 195, 198. On page 198 are a number of references to the 
use of the theory of elementary divisors in the study of differential equations. 
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The determinant of these n solutions can be written out in such a way 
that the elements above the principal diagonal are all zero. Thus the 
value of the determinant is : 


and the n solutions are therefore linearly independent. 

On account of the relation (10), each solution (11) of the system (8) 
will determine a solution of the original system (3). Accordingly we 
have n solutions of (3) which are linearly independent, for their de- 
terminant at any point is equal to the determinant of the A’s in (10) 
multiplied into the determinant of the solutions (11) for the same point, 
and neither of these determinants is zero. 

Suppose now we consider any multiple root of the characteristic de- 
terminant (5); for simplicity let us take 7,, and suppose that 


(12) 
80 that 
(18) 


is the multiplicity of the root r,. Then from (10) and (11) we see that, 
corresponding to this root, there are & solutions of (3) not involving log x 
and linearly independent : 


(14) yi? = Ay 


Therefore the constants A; , c, are linearly independent solutions of the 
equations (4) when r = 7,, as we readily see by putting the values of 


y;’** in (3). Now the equations (4), in this case, have only & linearly 
independent solutions, since there are only & elementary divisors corre- 
sponding to r,, and therefore at least one kth minor in (5) is not zero 


whenr=r,. If, then, we have any other solution of (3) of the form : 


= 0,2" (=1,2,...2), 


the constants C, must be linear combinations of the & sets of constants 
A; .,e We have then the following result : 


<=m 
$= 1,2,...8 


DUNKEL. — LINEAR DIFFERENTIAL EQUATIONS, 347 

The system of equations (3) admits n linearly independent solutions, 

such that, corresponding to each elementary divisor (r —1,)* of the 
characteristic determinant, there are e, solutions : 


l=e, 1 = 1 2 
Tk $ . 
(15) =2 3 A, .,1 (log 2) = 1, 2, 


If r, is a multiple root of the characteristic equation which furnishes s 
elementary divisors with the exponents e,, e, then the 


constants : 
, are s linearly independent solutions of the equations (4) when r = r,. 


SOLUTION OF THE CANONICAL SysSTEM IN THE GENERAL CasE 
BY Successive APPROXIMATIONS. 


We shall now return to the system of equations (1); and here again 
we shall make use of the linear transformation (10) to reduce the system 


to the canonical form: 
d 1 r 


The coefficients 4; are linear functions with constant coefficients 
of the coefficients a, , in (1). 

We shall now make use of the method of successive approximations to 
develop solutions of (16) about the point x = 0. It will be convenient 
to write the equations (16) in the form: 


d 1 ‘ 


The first approximation will be indicated by a third subscript 0, and is 
obtained as a solution of the system of equations resulting from (17) by { 
making the right side zero: 


* This reduction is used by Sauvage in the case of a system of equations with 
analytic coefficients. L. c., pp. 89, 90. 


§ 2. 
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d 1 r 


We have seen how to solve this system. Any one of the solutions 
that we have obtained may be used as the first approximation. This 
approximation having been chosen, we insert it in the right side of (17) 
and obtain the following relations for the first correction : 


d 1 


The right side is now a known function of x; and we have, conse- 
quently, a system of non-homogeneous linear differential equations to 
solve for z,,1- Having determined this first correction, it is inserted in 
the right side of (17), and the resulting equations are solved for the 
second correction. This process is repeated again and again, the relation 
connecting the gth and the (q + 1)th correction being: 


d 1 as 


Each equation (19) may be written: 


d 
whence : 
(20) 
= dx + fe S Sy as |. 


j=l 


Now writing out the value of z, ,1 4: in the same ways and substitut- 
ing it in the first integral of (20), we have: 


z z 


Ce, 
i=1 j=l 
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Now substitute in this result the value of 2, , 24:3 in the result thus 
obtained the value of z,,-5 .4:, etc. After a certain number of substitu- 


tions we have: 


Ce, 
Zz 


i=m j=e, 
+ flee. 41-+ g as | 


where Z < 1. 
When Z = 0 we have: 


2 2 


The lower limits ¢,, will be determined later to satisfy several con- 
ditions. 

We — choose at pleasure any one of the m elementary divisors, say 
(r —r,)**, and then select any one of the corresponding e, solutions of 
(18) S the first approximation. We shall take then for the first ap- 
proximation z, 9, the values given in (11) for a particular A; the inte- — 
gers x and ) will remain fixed for the solution we are now developing. 

For the development of the solutions corresponding to (r — r,)**, we 
shall make the following further assumption as to the coefficients ;7, 
Let us examine all the exponents of the elementary divisors (r — r,)°, 
which are such that Rr, = Rr,, where Rr, means “real part of r,, » 
and pick out one exponent, say e,, that is as great as any one in this 


special set of exponents; i.e., 
(23) 2 where Rr, = Rr,. 


The assumption is that the integrals : 


k= 1,3,...™ 
fi 


converge. If in particular r, is a simple root, and no multiple root has 
the same real part as it, then ey — 1=0, and this further restriction 
drops out. Or it might happen that r, is a multiple root, but that all the 
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exponents of the set (23) are unity, and in this case the restriction drops 


out also. 
The lower limits of integration in (22) will be determiued as follows : 


I. Rr, > Rr, 

(%,.=0 is L= 
Il. Rr,= Rr, 


HI. < Bre — 


where ¢ is a constant not zero, which will be determined more closely 
later (cf. p. 358). It will be proved in § 3, that even in the cases in 
which the lower limit is zero, the integrals converge. 

When all these conditions have been satisfied, we build hee the first 
approximation and the successive corrections the n infinite series : 


— k= 1, 2 m 


which will be proved in § 4 to converge and to form a solution of the 
system of equations (16). 

It will be convenient to consider in place of the functions z, ,, certain 
new functions ¢,,,, which will be defined by the following formulae : 


or Rr,=Rr,, k+xandl< 
or k=xand/<), 


(26) = (log Rr,= Rr, ktnxandL <l<e, 
x"* (log $x, A<lSe&. 


when Rr, + 


= dare 


For the case of g = 0, i.e., the first approximation, the ¢’s are certain 
_constants (cf. (11)); for all other values of g the formulae (26) define 
them as continuous functions of x so long as x is not zero: and we shall 
see later (Proof of Convergence, p. 358 et seq.) that each one approaches 
zero when x approaches zero. We shall therefore define each ¢,,,, 
when g > 0, as zero for x =0; and with this definition they will be 
continuous functions of zx in the whole interval 06. 

These ¢’s can be computed from the following recurrent formulae, 
which are easily obtained from (21), (22), (24), and (26) : 
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I. Rr> Rr. 


(27) 


De, 1, ¢+1 ee 
t=1¢ 


c = 

t integrations 


where h,;, is a positive integer whose value it will not be necessary to 
write out, noting, however, that : 


(28) hg S 
Il. Rr, = Rry. 
a) ISL. 
(29) = * = 
J * 
0 0 i=1 j=l 
t integrations 
b+ 
1 i=m je; 
i=1 j=l 
¢ 
ec) kaw, 15% 
tol 


; 
tintegrations 
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d) k=x, A<ISL. 


(32) tum = 


Ti=m j=e; 


i 
e) L<elS«e. 


(33) 
x, = dx, 2,941 (log dx 


Ill. Rr,< Rr,. 


Here we have the formula (29) again. 
In § 4 we shall consider the series: 


(34) > 


which are such that if we multiply each by its proper factor x re ( log)" aye, 
where h,, is given in (26), we obtain the m series (25). We tions 
reduce the proof of convergence of (25) to the question of the con- 
vergence of (34), and this last question will be settled by reference to 
certain formulae to be established in the next section. 


§ 3. 


Lemmas Concerninc INTEGRALS. 


We now prove a number of lemmas, which will be useful in the proof 
of convergence, and which also verify the statement that we have made 
that the integrals of the last section, in which the lower limit is zero, 
converge. 

Lemma I. Jf 6 ts a function of x, continuous in the interval 
O0<aSc, and such that |b||logx|', (¢ an integer > 1) ts in- 


| con 
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tegrable up to x= 0, then b can be integrated times from x = 0 as 
follows : 

t integrations 


Lemma II. Jf the conditions of Lemma I. hold, then: 


0 0 0 


integrations 


We will prove these lemmas by mathematical induction. They are 
true when ¢=1. Let us assume that they are true for a particular 


value of ¢, say ¢ = . 
Let X be any particular value of x in the interval 0< 2 <1, and — 
choose ¢ at pleasure such that 0 < «< X. Then we have: 


€ z € 
x 
— flog 2 [log ae 
e 
€ 
[og X1 ae 
e x 
S 10g 2 de + 10g 2 de 
: 


¢ x 
S ax + [log de = | B| [log 2 de, 
0 


VOL, XXXVIII. —23 
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Therefore, when ¢« approaches zero, the integral on the left side of all 
these inequalities converges, and we have: 


x x 
0 


Thus the two lemmas are proved when x < 1; and it is easily seen how 
to conclude the proof of I. in case ¢ > 1. 

Leva III. Jf 6 is a continuous function of x in the interval 
0<2x< ce, and its absolute value is integrable up to x = 0, and if: 


a1 
F(a) = b dx 
0 0 
t integrations 
where r is real and greater than zero, then: 


(87) flee. 
0 


When ¢ = 1 we have: 


z z z 


and so in this case III. is true. Assume that it is true for ¢ = ¢,; then 
it is also true for ¢ = ¢, + 1. For: 


< fewe][ fire] ne 


Therefore III. is true for all values of ¢. 


. 

| 

| 

| 

| 
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Lemma IV. Jf bis a continuous function of x inthe interval 0 <x <e, 


and if: 


t 
where r is real and greater than zero, then : 


(38) file! 
When ¢ = 1 we have: ' 


Assume that IV. is true for ¢=¢,. Then: 


<5 


Therefore IV. is true for ¢ = ¢, + 1, and the lemma is proved. 


Lemma V. Jf}, 1, and g, are defined as in IV. and it is further as- 
sumed that the absolute value of bi. ts lity up to x = 0, then: 


(39) limit 2” g, = 
r=0 
To prove this, let us choose a constant s such thatO <<s<r. Then: 


z 2 
ce ce 


= (r 
ce 


a~*|b| |de|, 


as we see from IV. by replacing r by r — s and 5 by ab, 
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If ¢ is a positive number chosen arbitrarily small, we can choose 7 so 


near 0 that: 
f 


and then X S 7 80 that: : 


x 


Then: 


sx feels fini <2< 


Therefore : 
|2" (x)| (O< x 4 X) 
and V. is proved. 
Lemma VI. Jf B is a function of x continuous in the interval 
and 
(40) limit B = 0; 


2 


(41) limit = 0. 


When ¢ = 0 it is obvious that (41) is true, for then: 
limit G(x) = limit B = 0. 
Assume now that (41) is true when ¢ = ¢,; then it will also be true 


for t=t,+ 1. For if ¢ is a positive number chosen arbitrarily small, 
we can choose y so near 0 that: 
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6, (@)|$$@+1 <1) 


(2) = f 1G, (2) de = fi 1G, de + fi G,, (x) de 
Hoge 


x 
(log — (log 


+ 
= § 1 — (P82) < 


Now choose X in the interval 0 < x < 7 80 that: 


1 1 


<$(4 +1) 


(log (#) 


1 
= Gr, () ae + Gt, (x) dx} S «, 


1 
og (2) = 0- 
Therefore VI. is true for all values of ¢. 
Lemma VII. Jf B ts a continuous function of x in the interval 


0< x Se <1, and such that it is not greater in absolute value than the 
constant N, and tf G, ts defined as in VI, then: 


© integrations 


fie. wf Lies... 


| | 
| 

and: 

| i 
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§ 4. 


Proor oF CONVERGENCE. 
The convergence of the series (34) will now be proved; and iu this 
proof it will also be shown that: 


(43) q>0. 
We shall use in the proof the following functions. A function of z, 
B, is chosen having the same properties of integrability as ;7, and such 
that: 
— k=1,2,...m 
44 B2 ). 
Next consider all the differences Rr, — Rr, which are not zero, and 
choose a positive constaut d such that : ; 
0<d<|Rr,— 


Then C shall be a positive constant such that : 


| 
(45) C2, C23 
We also define: a 
(46) M(x) = C B\log dex. 
0 


The point ¢ 1s chosen so that : 
(47) 0<ceSb, cS (where loge=1), M(c) < 1. 


This is the final determination of ¢ to which we have referred on page 
350; and this point ¢ will, from now on, mark one end of the interval 
for x. Instead of M(c), we shall, for the sake of brevity, write 
simply 

The convergence of the series (34) will be proved by showing by 
mathematical induction that the following inequalities hold for all 
values of 


k=m l=e, i=m jre, 
* We might, for instance,takeB==> 


47 
k=l l=1 i=1 j=l 
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b= 1,3..." 
(48) M* i= 1,2, . 


For g = 0 they are obviously true. Assuming them true for a special 
value of g, g = 91, we will consider the cases outlined in (24) in turn. 
I. From (27) we have: 


(49) 
t=l 
t integrations 


after replacing 
respectively by the greater values 
M", |log B 
(Cf. (28), (30), (32)). 


z Zz 
< M% | de| < M* OB|log | de| 


<M" M(2) < Me, 
From Lemma V. it follows that the limit of the right side of (49) is 
zero, when x approaches zero; for we have assumed that B| log x |°«—* 


is integrable up to x = 0 (cf. page 358), and therefore B| log x |**™ 
must also be integrable up tox =0. So (43) is verified for g = g, + 1. 
II. From (29) we have for case (a) : 


t—l z 
0 0 


t integrations 


— 
<mu*> f B|log dee (Lemma II.) 


t=l 
aig . 


| 
a 
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since the greatest value of e, — A +¢— 1 is 


— Mae 
G1) MOL J Blog de = M2) 


< M"M (2) < 


It is obvious that the above result holds also for II. (c) (31). From 
(32) we shall obtain in the same way for II. (d) : 


0 


Mn" 


| atl jlog <M M% M (a) < Met 


In all three sub-cases (a), (c), and (d), it is easily seen that (43) is 
true for g=qi+1. We have now left of case II. the sub-cases (6) 
and (e). 

From (30) we have for case (6) : 


(53) 


1 
Lf? |dx|.. of fare ia 


L) integrations 


e € c 

t integrations 


In the first part of the bracket we have replaced 1u,z,¢,41! by 
une (x), using inequality (51). 
We will consider the two parts of (53) separately. 
(54) z 
(1 — L) integrations 
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= 1 1 


{Integrations 
L 
x f B|log x ae | (Lemma VII.) 
t=1 
t=I—L 


Therefore : 
1 +1 l 1 
I de, |< < af = M% <M a+ 


The limit of the left side of (54) for x= 0 is zero by Lemma VI, 
while the same thing is true of the left side of (55) from the inequali- 
ties. Therefore (43) is true in this case for g = q, + 1. 

From (33) we have for case (e): 


(56) 
1 
e e 


(l—L) 


t 


In the first part of the bracket we have used the inequality obtained 
from (52) : 


L 
[log x 6. S (2). 


Now the only difference between the inequalities (53) and (66) ie is in 
the power of |log a| outside the brackets; and, since | log x| git 
| log x |“, all the results that we have obtained from (53) will follow 
also from (56). “Then for all the sub-cases of II., (43) and (48) are 


true for g = q, + 1. 


» 
‘ | 
a | de | 
il 
| 
i 
{ 
| 
{ 
| 
| 
q 
| 
| 
4 
a 
| 
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III. Here we have: 


1, 


t=l z z 
Sf Lae... f BM ae 
t integrations 


< M% f Blog x dex (Lemma III.) 


t=1 R (Te r,) 


< \log de < 
0 


It is obvious that (43) is true in this case for g = g, + 1. 

This completes the proof of the inequalities (48). 

From (48) we now infer at once that ¢,,,, (¢ >) approaches zero at 
the point x = 0. 

The inequalities (48) furnish the sufficient condition of Weierstrass 
for the absolute and uniform convergence of the series (34) in the in- 
terval OS 2<ec. Accordingly the series (25) converges absolutely 
and uniformly in any sub-intervalO Cel 

It remains now to show that we have a solution in the system 
of functions z,,. For this purpose let us select any sub-interval 
Since the point = 0 is excluded, the coefficients 5; 
are continuous in this interval; and in it the system of differential 
equations (16) is satisfied by the functions z,, For if we multiply the 
series (25) for z,, by 5;”, we obtain the absolutely and uniformly con- 
vergent series : 


q=n 
> q=0 


Taking the sum of such series for all values of ¢ and j we have: 


i=m j=e, 


i=m 
(57) > = > > > 
i=1 j=l 


i=l 


q—0 j=l 
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where the last transformation is valid since we are dealing with abso- 
lutely convergent series. Also: 


(58) 


(59) 


Adding (57), (58), (59), and changing slightly the summation on the 
right, we have : 


1 
1 


1 


i=1 j=l 


The series on the right is an absolutely and uniformly convergent 
series of continuous functions. If we replace the terms by their values 
given in (18) and (19), we can also write the series in the form: 


d = d 
dx dx dx 


which is the series of derivatives of the terms of (25). From this it 
follows that, if we differentiate the series for z,, term by term, we shall 
obtain an absolutely and uniformly convergent series of continuous 
_functions ; and therefore z,, has a continuous derivative at each point 
of our sub-interval, which is precisely (61), and this, as we have seen, 
is the same as the right side of (60). We can therefore write (60) in 
the form 


and now giving & and / all possible values we have precisely the system 
of equations (16). We have, then, the following result: 


In any sub-interval of 0c, not including the point x = 0, the functions 
z,,, represented by the series (25) are continuous in x, have continuous 


|| 
1 1 
q=0 
q=0 
| 
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first derivatives with regard to x, and satisfy the canonical system of 
equations at every point. The terms of the series for z,,, are given by the 
formulae (26), in which the functions $,,, are continuous in x through- 
out the whole of the interval 0b, and reduce to zero when x = 0, except 
wheng=0. The functions ¢$,,,,9 are constants, zero in all cases but the 
following : 
1 


§ 5. 


Linear INDEPENDENCE OF THE SOLUTIONS OF THE CANONICAL 
System. 


We have shown that corresponding to each elementary divisor 
(r — rx)°*, there are e, solutions of (16) obtained by giving A the values 
1, 2,...¢ 3; and for the development of these solutions we have re- 
quired that | | log ax shall be integrable up to x = 0, where 
ex is determined for the root 7, by the condition (23). We have, then, 
n solutions which may be written as follows: 


or Rr, = 

or k=«xand/<), 
(62) = (log Rn=Rr, kt L 
= (log x) A<lSe. 


Rr,t Rr, 


where the functions are continuous in x and 


1 
(I—A)! ASIse, 


(«= 1, 2,...m) (A = 1, 
It is worth while to note three facts in regard to the z’s, which will be 
useful later on : 
I. 2 does not involve log z explicitly whenever A > 7; 


(64) IL, init = 0 when Rr, > Rr,; 


(63). 


Ill. ini (log x) = 0 in all cases except the 


1 


“one, and I= and then the is 


| | 

| 
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We will now show that the z solutions (62) are linearly independent. 
Suppose they were not and that there were n constants C,,, not all zero 
such that: 


x=m A=e 
k=1,2,...m 


It will be convenient to suppose that our notation is such that: 
(66) Rees... Rte. 


Consider first those equations of (65) for which Rr, = Rr, We 


have: 


(i= 1, 2, 


Now let /= 1 in (67), and consider the limit of each term for any 
given value, within the range indicated, for &. For each term in which 
Rr, > Rr, = Rr, the limit is zero by II. We have left, now, only 
the terms: 


(68) Rr, = Rr, = Rn, 


According to I. no logarithms appear explicitly in (68), and we can 
write: 


Now by (63) the limit of all such terms is zero except in the one case 
«x == k and A = 1, and for this term we have: 


So in the case of 7 = 1 the limit (67) turns out to be C,, when 
we evaluate the limit term by term. Now this is impossible unless 
C,,1 = 0; and so we must have, writing now « in place of k: 


(71) C.,1 = 0 Rr, = 
Now consider in the same way the cases of (67) in which 7 = 2; and 


choose any one of the values of / indicated. Here again by II. the limit 
of éach term is zero when Rr, > Rr,; and we have left the terms: 


(72) Rr, = Rr, = Rr, 
r 2 2. 


| 
q 
f 
if 
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Now by I. no logarithms appear explicitly in the terms (72), for the 
only case in which they could occur would be for A = 1, but by (71) 
such terms do not appear. We can then write each term of (72) 


and by (63) the limit of each is zero except in the one case: 


(74) limit = limit = 


We have then, reasoning as before, and including the previous result, 


2=0 Rr. = Rx. 


(75) = C, 
Now the same reasoning can be applied until we have finally 


(76) Cx, = 0; Rr, = Rr,, 
A>1 


and there are left in (65) only those solutions for which Rr, > Rr. 
Having thus disposed of all the 7’s in (66) such that Rr, = Rr,, we 
consider the next set of r’s whose real parts are all equal and as small as 
any other in the new set of r’s, and show, in exactly the same way, that 
the corresponding constants, C,,, are zero. Continuing in this way, we 
finally reach the result that all the constants in (65) are zero, and that 
the assumption made as to the dependence of the n solutions leads to a 
contradiction. Zhe n solutions are therefore linearly independent. 


§ 6 


Return THE CANONICAL SysTeM TO THE System. 


We will now return to the original system of equations (1). Each 
solution of the canonical system (16) determines a solution of the original 
system (1) on account of the relation (10). We have then the following 
n solutions of (1) : 


k=m | 


(77) > D> 


i=l ove Cx 


Now these n solutions are linearly independent, for the z solutions are 
linearly independent, and the determinant of the constants A,,, is not 


zero, 
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Now consider the following limit : 


(78) limit (log a) 
k=m I=e, 
=limit > > (log a) 
k=l 
On the right, if we take the limit of each term in the summation, we 
shall find that it is zero in every case except the one in which k = « and 
l= e, (cf. (64) III.). We have, then: 


1 


So we can write as a set of solutions of (1), corresponding to the ele- 
mentary divisors (r — 7.) : 

(80) = (log x) 1, pee 
where the functions vr are continuous in the interval 0 < x < ¢, and 
such that : 


The constants A,;,, are determined independently of the functions 
a, in (1); and therefore all that we have said on page 346 in regard 
to certain sets of them as linearly independent solutions of the equations 
(4) in the special case of a, ; = 0 holds equally well here. 

In order to obtain solutions corresponding to (r -- r,)°*,* we assumed 
in the treatment of the canonical system that |8;4||log 2|°*—* was inte- 
grable up tox=0. Now, since the coefficients bj’, are linear functions 
with constant coefficients of the coefficients a,,, it will be sufficient, in 
order to obtain solutions of (1) corresponding to (r — 7,)°*, to assume 
that |, ,| |log a |x? is integrable up to x = 0 for all values of ¢ and 7. 

Our results may be stated as follows: Jf (r —1,)“* ts an elementary 
divisor of A(r), and if we consider all the elementary divisors (r —1r,)"* of 
A(r) such that Rr,= Rr,, and denote by ex that exponent which ts as 
great as any other exponent in this set of elementary divisors, and assume 
that |a;,| |log «|°*— is integrable up to x = 0; we can develop e, solutions 
of (1): 


* Cf. p. 849. 


\ 
| | 
| 
| 
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x, « = 1 2 N 


where 5" is continuous in the neighborhood of x = 0 and 
1 


If the root r, furnishes s elementary divisors : 
then the constants : 


Ai,x,e,5 = 1, 2 eee n) 


are s linearly independent solutions of the equations (4) when r = rx. 
The n solutions of the differential equations (1) that we obtain when 


«= 1, 2,...m are linearly independent. 


§ 7. 


Tue Homoceneous Linear DIFFERENTIAL EQUATION OF THE 
nth ORDER. 


We shall consider homogeneous linear differential equations which can 
be written as follows : 


(82) + ( saa tees 


dx" 
% 1 


in which 1, are constants, and p,, pe,...p, are functions 
of the real independent variable x, continuous in the interval 0 < 2 < 4, 
and such that their absolute values are integrable up to x = 0; in short, 
these p’s are to have the same properties as the functions a, , in (1). 

This equation can be replaced by a system of linear differential equa- 
tions by the following substitutions : 


(¢=0,1,...n—1). 


I 


(83) 
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We thus get the system of differential equations: 


1 


x 
(84) 


The characteristic determinant of (84) is : 


r—n+1+p, 
—1 


(85) 
0 
0 


The adjoints of the first line are: 
(86) r(r—1)...(r—n+ 2); 
r(r—1);3 7; 1. 
The characteristic equation is, then : 
(87) 4... 
(7 — 1) + Pua? + = 0.* 
There is always one first minor of (85) which is not zero, the adjoint 
of »,; and therefore if r, is a multiple root of (87), the only elementary 
divisor of (85) corresponding to r, is (r — 7,)“*, where e, is the multi- 
plicity of the root r,. 
We have seen that, corresponding to the elementary divisor (r — r,)**, 
there are e, linearly independent solutions of (84) of the form (80). 
Here the constants Aj,x,¢, must satisfy the set of equations : 


(88) (te —M +1 + + + An ce, = 0 
= 2, 8). 
This system of equations has essentially only one solution, namely : 
(89) Aine, = (¢=1,2,...n—1), 


An, — Pr» 


* (87) is also called the indicial equation of (82) for the point x = 0. 
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where p is a constant not zero. We may divide each solution by the 
jw?’ and the resulting solutions may be 


(e 
written : 


A ¢mi,3,...8 


where =1.e(te—1) ++. +14 4), 
Ky A 
= 1. 


By means of (83) we can now return to solutions of the equation (82) 
with the following result : 

If r, is a root of the characteristic equaticn, and | p;| | log lis 
integrable up to x = 0, where ex > ¢, for all k’s such that Rr, = Rr, 
e, being the multiplicity of the root r,, then the equation (82) has e& 
linearly independent solutions which may be written with their first n — 1 
derivatives : 


corresponding number 


a (log 2)" 


(91) 
“da 


Aw 1,2,...% 


where the functions Ee”. are continuous in the neighborhood of x = 0 
and, 


(92) 


yA 
= 


= fa (te — 1) — + 2). 


Even for the equation of the second order this theorem does not 
merely give the results of Professor Bécher’s paper above quoted, but 
goes a step farther, since in the case in which the two roots of the 
characteristic equation are equal, we require merely that: 


b b 


converge, while Professor Bécher’s method made it necessary for him to 
require that 


b b 
2) fog 2)*| ps| ae 
0 0 


converge.* 


* L.c. p. 48. The function in 


